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ABSTRACT 


A complete  approximate  set  of  equations  is  developed  to  describe  fire  spread 
over  a surface  and  its  resultant  energy  release.  Wall,  floor,  and  ceiling 
orientations  are  considered.  The  needed  model  data  are  couched  in  terms  of 
available  test  method  results,  e.g..  Cone  Calorimeter  and  LIFT  apparatuses. 
Several  applications  are  presented  to  show  how  energy  release  rates  can  be 
predicted  and  how  well  they  represent  real  data  from  full-scale  and  model  room 
lining  experiments. 


Key  Words:  Heat  release,  flammability  tests,  fire  prediction,  linings,  room 
fires . 


1 . INTRODUCTION 

Traditionally,  fire  test  methods  have  been  used  to  provide  relative  rankings 
of  material  and  product  performance  from  which  officials  have  defined  safety 
requirements.  This  process  has  been  empirical  relying  on  judgement  and  fire 
experience.  Although,  one  would  expect  fire  safety  to  be  consistent,  each 
regulatory  agency  or  government  has  laid  out  its  own  tests  and  procedures 
creating  an  array  of  different  results.  Test  methods  stand  more  on  their 
longevity,  than  on  their  scientific  merits.  This  is  an  unsound  technical 
situation  which  can  compromise  safety.  The  solution  to  this  state  of  affairs 
is,  however,  not  readily  accessible  since  a technical  consensus  has  not  yet 
been  developed  from  the  science  of  fire.  Hence,  a complete  scientific 
alternative  is  not  yet  possible,  but  a framework  for  an  alternative  can  be 
developed  in  terms  of  available  scientific  information  and  methods.  We  will 
present  such  a framework  and  demonstrate  its  utility  and  level  of  accuracy. 
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Theories  of  flame  spread  and  new  widely  used  test  apparatuses  capable  of 
providing  engineering  fire  property  data  make  it  feasible  to  construct  a 
framework  for  fire  hazard  assessment.  The  test  apparatuses  include  the  Cone 
Calorimeter  [1],  and  the  Lateral  Ignition  and  Flame  Spread  Test  (LIFT)  [2], 
both  of  which  have  been  recently  established  as  ASTM  standard  methods  and  are 
being  considered  in  ISO.  We  shall  not  review  the  details  of  these  apparatuses 
here,  it  only  suffices  that  they  can  produce  the  needed  data  in  the  analysis 
to  follow.  That  analysis  will  be  confined  to  fire  growth  on  flat  materials  in 
applications  of  floors,  walls  and  ceilings.  We  shall  formulate  a model  to 
compute  the  "output"  of  the  fire  growth  as  a function  of  time  in  terms  of  the 
material  orientation  and  properties.  The  "output"  can  be  defined  as  ap- 
propriate for  the  hazard.  In  the  exercises  to  follow  we  shall  address  the 
output  as  energy  release  rate  representative  of  the  hazard  of  rapid  fire 
growth  or  perhaps  flashover. 

The  model  to  be  presented,  we  emphasize,  represents  only  a framework  for  a 
more  complete  and  precise  solution  to  the  appropriate  characterization  of 
materials  through  fire  test  methods.  But  once  formulated,  we  tested  the  model 
against  various  data  sets  for  full  scale  behavior  of  fire  growth  on  lining 
materials  in  rooms.  Our  results  and  our  model  should  serve  as  support  for 
other  analyses  which  have  used  similar  property  data  to  predict  the  flashover 
times  for  some  of  these  experiments  [3-5].  However,  we  believe  that  our  model 
is  less  empirical  and  more  physically  complete  than  these  previous  analyses 
[3-5],  and  the  accuracy  of  the  end  results  is  comparable.  Indeed,  it  suggests 
that  we  all  have  selected  the  key  properties  in  our  analyses  which  illustrates 
the  value  of  the  new  fire  property  tests  in  supplying  useful  and  versatile 
engineering  data.  Hence,  future  models  should  be  able  to  build  improvement  on 
this  base.  Finally,  our  model  will  be  capable  of  coping  with  changing 
environmental  conditions,  but  we  shall  not  implement  this  in  any  dynamic  sense 
by  coupling  its  fire  growth  to  the  changing  conditions  of  a room  involved  in 
fire.  However,  we  will  examine  some  effects  of  pre -heating  during  flame 
spread  that  could  come  as  a result  of  thermal  feedback  due  to  a confined  fire. 


2.  THEORETICAL  DEVELOPMENT 

Initially  the  purpose  of  the  analysis  to  follow  was  to  illustrate  how  fire 
property  test  data  could  be  used  in  modeling  to  yield  meaningful  quantitative 
results  for  assessing  fire  safety.  As  will  be  seen,  the  use  of  the  analysis 
is  capable  of  producing  very  accurate  results . This  probably  resulted  from 
our  approach  of  keeping  the  model  simple  but  complete.  We  have  included  the 
essential  elements  of  fire  growth;  namely,  ignition,  energy  release  rate, 
burning  time,  and  all  relevant  modes  of  flame  spread  for  the  surface  orienta- 
tion considered.  Also,  we  have  included  effects  of  the  ignition  fire  source, 
and  have  considered  parametrically,  the  effects  of  thermal  feedback.  Detailed 
transient  characteristics  and  room  smoke  layering  have  been  ignored,  but  these 
effects  could  easily  be  included  in  a more  refined  analysis  and  their 
contributions  assessed.  However,  it  can  be  argued  that  their  impact  is  likely 
to  be  secondary,  but  ultimately  needs  to  be  examined. 
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Fire  Properties: 


The  term  "fire  properties"  has  been  used,  but  we  need  to  define  it  in  context. 
Generally,  engineering  property  data  constitute  physical  and  chemical 
parameters,  usually  arising  from  thermodynamic  principles  or  rate  process 
models.  Although  we  could  couch  our  analysis  in  terms  of  such  true  proper- 
ties, we  choose  to  use  more  directly  measurable  quantities  since  fire  modeling 
is  not  yet  complete  enough  to  deal  with  complex  materials.  Hence,  our  "fire 
properties"  will  not  necessarily  be  valid  engineering  properties  and  will 
depend  on  environmental  conditions.  We  will  identify  these  "fire  properties" 
as  we  develop  our  theory.  As  an  introduction,  we  list  them  now: 

1.  Q" , average  peak  heat  release  rate  per  unit  area 

This  value  is  derived  from  a device  like  the  Cone  Calorimeter  [1]  for  which 
the  transient  data  have  been  averaged,  and  evaluated  at  a specified  irradiance 
level.  It  is  illustrated  in  Figure  1 in  which  the  area  under  the  actual  curve 
is  approximately  equal  to  the  area  under  the  rectangle.  We  have  not  used  a 
standardized  prescription  for  determining  the  average  energy  release  rate  and 
burn  time.  This  was  done  in  a consistent  manner  by  approximating  the 
significant  burning  region  of  the  curve  by  a square  wave.  We  selected  an 
average  maximum  for  Q"  and  then  estimated  t^^  to  maintain  the  total  energy. 

2 . tjj  , burn  time 

This  value  is  the  width  of  the  approximating  rectangle  in  Figure  1.  This  also 
corresponds  to  a given  heat  flux,  and  is  derivable  from  the  Cone  Calorimeter. 

3.  ignition  time 

This  value  is  the  time  for  piloted  ignition  at  a given  heat  flux,  specifically 
radiative.  The  heat  flux  level  should  be  selected  to  correspond  to  the  flame 
heat  flux  level  appropriate  to  wind-aided  flame  spread.  This  will  be  defined 
more  specifically  later.  This  "property"  can  be  derived  from  a number  of 
ignition  apparatuses  without  any  loss  in  uniqueness  since  it  should  not  be 
very  dependent  on  the  apparatus  (e.g..  Cone  Calorimeter  [1],  [2])  provided  the 
radiant  emission  spectrum  of  the  apparatus  is  similar. 

4.  flame  heating  parameter  in  opposed  flow  spread 

This  "property"  is  directly  derivable  from  the  LIFT  [2,6]  standard  procedure. 
It  represents  the  flame  heat  transfer  and  opposed  flow  velocity  effects.  In 
the  LIFT  mode,  the  flow  velocity  corresponds  to  a natural  convection  condi- 
tion, and  experience  suggests  that  the  $ is  not  strongly  sensitive  to  surface 
orientation. 

5.  Tg  min’  minimum  temperature  for  flame  spread 

In  opposed  flow  flame  spread  it  is  found  that,  in  general,  materials  require  a 
minimum  heating  level  to  achieve  a surface  temperature  required  for  spread. 
Below  this  temperature  spread  is  not  possible.  Above  this  temperature  and  up 
to  the  ignition  temperature  (T^^)  of  the  material  is  the  domain  of  surface 
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temperatures  where  surface  flame  spread  theory  applies.  For  external  heating 
conditions  that  lead  to  the  attainment  of  the  ignition  temperature,  a gas- 
phase  flammable  mixture  can  result  with  flame  propagation  still  possible  as  is 
a "gas -phase"  phenomenon  in  contrast  to  surface  spread  that  we  are  considering 
here.  Although  this  property,  can  be  derived  from  the  LIFT  for 

opposed  flow  spread,  it  is  also  required  but  not  available  for,  wind-aided  or 
upward  flame  spread.  For  example,  it  appears  that  a pre-heated  surface 
temperature  in  excess  of  80“ C is  necessary  to  cause  sustained  upward  flame 
spread  on  wood  particle  board  [7] . 

Equations : 

In  developing  the  governing  equations  we  shall  consider  all  materials  to  be 
thermally  thick,  i.e.  no  effect  of  thickness,  and  laminated  materials  and 
substrates  are  treated  as  homogeneous,  reflective  of  their  bulk  properties. 
The  same  corresponding  equations  will  apply,  without  distinction,  to  horizon- 
tal (floor) , lateral  and  downward  wall  spread  as  opposed  flow  spread;  and  to 
upward  wall  and  ceiling  spread  as  wind- aided  spread.  One  might  argue  with 
this  lack  of  distinction  but  the  opposed  flow  case  for  these  different 
orientations  has  some  support  [8]  whereas  the  extension  of  the  upward  spread 
theory  to  the  ceiling  is  only  a convenience  for  now  which  needs  further  study. 
For  all  orientations,  as  illustrated  in  Figure  2,  symmetry  is  invoked,  and  the 
wind- aided  coordinate  is  designated  by  y and  the  opposed  flow  coordinate  by  x. 
Spread  on  the  floor  (Fig.  2a),  and  ceiling  (Fig.  2b)  is  radially  symmetric, 
and  on  the  wall  (Fig.  2c)  is  approximated  as  a region  developed  from  two 
linear  (x  and  y)  spreads  away  from  an  initially  ignited  rectilinear  zone  of 
area  2'x.^y^  . The  region  between  the  extended  pyrolysis  wall  fronts  (Xp  , yp ) , 
and  the  burn-out  fronts  (Xjj,yjj)  is  formed  by  straight  line  approximations  as 
illustrated  in  Figure  2c. 

The  objective  of  our  analysis  is  to  compute  the  energy  release  rate,  Q(t) , as 
a function  of  time  for  the  orientation  shown  in  Figure  2 or  their  combination. 
For  energy  release  we  write  the  following  approximation  based  on  Figure  1. 

Q(t)  - JJ  Q"  dx  dy  = q"  Ap  (1) 

where  Ap  is  the  pyrolysis  or  volatilizing  area. 

For  the  floor  or  ceiling  cases  it  can  easily  be  shown  that 

Ap  - X (Xp2  - Xp2)  or  IT  (yp2  . (2) 

For  the  wall  case, 

= 2 [ypXp  + (Xp-Xp)yp  + 1/2  (yp -y, ) (Xp -x^  ) ] for  t < (3a) 

or 

'S,  “ 2 [ypXp  + Xpyp  + 1/2  (yp -yp  ) (Xj, -Xp  ) ] for  t > t;,  (3b) 

- 2 [yfcXp  + Xi,yp  + 1/2  (yi, -yp  ) (Xt -Xp  ) ] 

In  order  to  formulate  the  rest  of  the  analysis  we  write  out  and  solve  the 
equations  for  the  pyrolysis  and  burnout  fronts. 
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Wind-aided  Solution: 


The  upward  or  wind-aided  solution  is  developed  [9]  from  the  equation 


<iyp  y£  - yp 


dt 


(4) 


where  is  the  position  of  the  flame  tip  with  - y defining  the  forward 
heat  transfer  region  of  the  flame  identified  with  heat  flux,  q^' . Equation  (4) 
is  a quasi-steady  model  which  uses  the  steady-state  flame  spread  result 
successively  over  time  and  the  forward  heat  transfer  region  changes. 
Mathematically,  this  solution  is  continuous.  This  flame  heat  flux  will  be 
taken  as  30  kW/m^  in  our  analysis  for  all  materials. 

The  ignition  time  is  given  as 

)r  kpc  (Tjg  - 
‘"is  4 (qj)2 

where  T^^  is  the  ignition  temperature,  kpc  is  the  thermal  inertia,  and  T^  is 
the  surface  temperature  of  the  region  y > y^ . As  can  be  seen,  the  "property” 
t^g  can  be  related  to  more  valid  engineering  materials  properties  (kpc  and 
T^g)  and  other  factors. 

Finally,  the  flame  length  can  be  estimated  by  a linear  approximation  to  wall 
fire  data  [9]  as 

yf  - yb  = ''f  [Qo  + Q (yp  - yb)l  (6) 

where  y^,  is  taken  as  identically  zero  for  t < t^^  , and  is  the  energy  release 
rate  per  unit  width  of  the  ignitor  flame  needed  to  start  the  process.  This 
ignitor  flame  is  assumed  to  provide  a uniform  heat  flux  q^  = 30  kW/m^  over  the 
region  which  is  taken  to  ignite  at  t = t^g  in  Eq.  (5),  or  correspond- 

ing to  t^g  in  Figure  1.  The  value  for  k^  is  approximately  0.01  m^/kW.  In 
Equation  (6),  the  term  is  set  to  zero  after  t = tj^  . This  is  done  after  t = 
tjj  since  the  initial  region  as  determined  by  the  flame  extent  due  to  is  no 
longer  burning,  and  subsequently  the  ignitor  flame  would  become  discontinuous 
from  the  material  flame  ahead  of  this  region.  This  ignitor  energy  is 
consistently  included  in  the  overall  heat  release. 


The  burnout  front  can  be  approximated  as 

dyb  yp  - yb 


dt 


for  t > U 


The  initial  conditions  are  given  as: 

t = 0 , yp  = y^ 

^ ^ . yp  - yp(tb)  . yb 


(7) 


(8a) 

(8b) 


and 
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(9) 


The  solutions  to  the  above  are  given  below: 

V = (— ) e - 1/a 

for  0 < T < t^/t^  where  r = t/t^  , r?  =»  Yp/Yo  . ^ ^ ^fQ"  ” ^ Yo  = ^fQo  * 


This  assumes  that  the  initial  pyrolysis  length,  y^ 
height  of  the  ignition  flame 


is  determined  by  the 


For  r > (tjj/t^  ) ^ ’’b  • Eqns.  (4)  and  (7)  are  subtracted  and  solved  to  yield 


Yb 


where  b = a - l/r^^ 


c e 


b(r  - r^) 


(10) 


, /a+lx  r 1 1 

and  c = (— ) [e  ® - 1] 

3. 


From  Eqn.  (7)  it  follows  that  i}^  -=  Yb/Yo  given  by 


^b  ^ 1 + 


br. 


. 1) 


(11) 


The  parameters  a and  b must  be  greater  than  zero  in  order  for  the  upward  fire 
spread  to  accelerate.  For  values  less  than  zero,  the  spread  will  eventually 
stop.  This  bimodal  behavior  is  a distinct  characteristic  of  upward  flame 
spread  which  is  conditional  on  the  form  of  Eqn.  (6)  . This  clearly  shows  the 
role  of  the  energy  release  rate  per  unit  area  and  its  duration  in  time,  i.e., 
Q"  and  t^  . 

Opposed  flow  solution: 

From  the  LIFT  procedure  [2,6]  the  governing  equation  for  the  pyrolysis 
position  is 


dt 


$ 


for  T > T 


kpc  (Tig  - T^)^  "s.min- 


(12) 


This  can  be  rewritten  as 


TT  ^ 


dx  4 (q^) 


dt 


'ig 


The  dimensionless  solution  | 

I = + 1 


Xp/Xg  follows  as 


(13) 


(14) 


where 


TT  X_ 

4 (q^ 
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The  burnout  solution  analogously  follows  from  Eqn.  (7),  i.e. 


and  it  can  be  shown  that  = ^b/^o  given  by 
lb  = ^ 1 ^ ’’b  • 

This  completes  the  solution.  These  equations  will  now  be  applied  to  specific 
applications.  For  each  case  we  shall  appropriately  apply  the  equations  by 
accounting  for  the  effects  of  geometry  and  thermal  feedback  as  a parameter. 
For  example,  a corner  wall  and  ceiling  fire  will  be  considered  as  an  extended 
vertical  wall  representative  of  the  actual  upward  and  wall  ceiling  jet  region 
plus  a quadrant  of  a circular  ceiling. 

To  review,  the  procedure  for  implementing  this  solution  is  as  follows: 

1.  Assemble  material  data  namely,  Q" , and  T^  ^in  • 

2.  Prescribe  the  ignition  source,  . 

3.  Since  no  room  fire  model  is  coupled  to  this  analysis,  one  must 
set  Tg  to  a value  representative  of  the  particular  room 
configuration . 

4.  Compute  the  pyrolysis  area  for  the  surface  configuration 
using  the  formulas  for  Xp  and  yp . 

5.  Compute  the  fire  product  output:„rate , e.g.  rate  of  any  energy 
release,  by  multiplying  Ap  by  Q ’ 


3.  APPLICATIONS 

Several  applications  of  the  analysis  will  be  presented  for  cases  where 
scenario  results  for  total  energy  release  exist  along  with  a complete  set  of 
fire  property  data. 

3.1  Full  Scale  Room  - Corner  Tests 

A series  of  full  scale  room  lining  fire  tests  were  performed  at  the  Swedish 
National  Testing  Institute,  Boras,  Sweden  [10].  The  test  scenario  consisted 
of  lining  the  walls  and  ceiling  of  the  room  (2.43  m by  3.66  m by  2.43  m in 
height  with  an  opening  2 m by  0.9  m)  with  a given  material.  A 0.17  m by  0.17 
m sand  burner  was  placed  in  the  corner  of  the  full  scale  room.  Initially  a 
100  kW  fire  was  allowed  to  burn  in  the  corner.  If  this  fire  did  not  cause 
flashover  in  10  minutes,  then  the  burner  output  was  increased  to  300  kW. 

Small  scale  fire  test  data  were  obtained  by  other  researchers  for  the 
materials  used  in  the  Boras  study  [5,11,12]  and  were  used  to  develop  Table  1. 
These  values  are  used  for  the  model  inputs.  In  the  calculations,  ceiling  fire 
spread  is  accounted  for  by  having  the  ceiling  ignite  when  the  pyrolysis  front 
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reaches  the  ceiling.  The  ceiling  spread  follows  the  wind-aided  spread  along 
the  wall  ceiling  interface.  Thus,  the  pyrolyzing  ceiling  area  is  given  as  the 
quadrant  of  the  ceiling  with  the  radius  equal  to  the  distance  from  the  corner 
to  the  wall  ceiling  interface  pyrolyzing  front.  Lateral  spread  was  deemed  to 
occur  only  when  the  wall  surface  temperature  is  at  the  material's  T^  or 

greater.  The  model  was  run  with  T^  at  25 *C  and  80" C.  While  at  the  start  of 
the  experiment  the  surface  temperature  is  uniform  (approximately  25 "C) 
throughout,  allowing  T^  to  be  an  input  parameter  shows  the  sensitivity  to 
feedback  and  the  effects  which  increased  lateral  spread  has  on  the  results. 
This  is  justified  because,  while  early  on  in  the  test  T^  of  25 "C  is  more 
representative  of  the  room  conditions , the  upper  layer  temperature  is  always 
increasing  and  radiative  feedback  is  raising  wall  and  ceiling  temperatures. 
We  selected  burner  flame  heights  (y^ ) of  1.3  m and  2.1  m for  the  100  and  300 
kW  settings . We  determined  y^^  values  from  an  axisymmetric  fire  plume 
correlation  realizing  that  corner  flames  would  be  taller  and  ignoring  ceiling 
interaction.  For  cases  that  stopped  propagating  before  10  minutes,  the 
analysis  was  restarted  at  t=600  s from  a new  origin  at  y^  (600  s)  . If  y^j(600 
s)  is  greater  than  2.1  m,  the  ignitor  flame  height,  we  arbitrarily  selected 
0.3  m as  a new  y^  in  order  to  initiate  spread  again.  Also  for  t>  600  s,  we 
selected  T^  = 250" C based  on  the  maximum  gas  temperatures  measured  [10]  . 
Model  predictions  and  experimental  results  for  the  time  to  a 1 MW  fire  are 
given  in  Table  2.  Figure  3 shows  the  experimental  rate  of  heat  release  as  a 
function  of  time  (solid  curves)  and  model  calculations  (dashed  curves)  with  T^ 
= 80"C  for  three  of  the  materials.  Appendix  A contains  plots  of  experimental 
rate  of  heat  release  and  model  calculations  for  all  13  materials.  The  results 
are  comparable  to  those  of  Karlsson  (who  only  considered  6 of  the  materials) 
[5],  but  appear  to  be  distinctly  better  than  the  results  of  Wickstrom  and 
Goransson  [3]. 

3.2  PMMA  Room  Wall  Test 

A PMMA  wall  fire  test  was  performed  at  NIST  [13].  The  wall  opposite  the 
doorway  of  an  enclosure  of  roughly  the  same  dimension  as  the  Boras  room  was 
lined  with  12.7  nun  thick  clear  PMMA,  and  a small  line  burner  0.3  m wide  was 
placed  at  floor  level  in  the  center  of  the  wall.  The  igniter  flame  height  was 
nominally  0.15  m.  The  pyrolysis  area  as  a function  of  time  was  recorded  from 
video  records  and  imbedded  thermocouples.  Initially,  the  fire  spreads  up  the 
wall,  and  across  the  wall-ceiling  interface  (ceiling  jet).  This  mode  of 
spread  is  characterized  as  wind- aided  spread.  As  this  is  occurring,  lateral 
spread  is  extending  the  initial  pyrolyzing  width  of  0.3  m.  Appendix  B shows 
the  measured  surface  temperatures  at  different  position  and  the  estimated 
pyrolysis  area. 

The  model  input  values  for  this  PMMA  wall  fire  are  burning  rates  of  6,  7,  and 
8 g/m^s  (multiplied  by  25  kJ/g  for  rate  of  heat  release),  160  s for  t.  (this 
assumes  a value  of  380"C  for  , 1.02  (kW/m^K)^s  for  kpc,  and  q^*  of  25 

kW/m^ ) , tjj  of  1000  s,  and  $ of  14.4  kW^/m^  . The  minimum  temperature  for 
lateral  spread  is  approximately  25 "C  for  this  material.  With  these  input 
values,  the  pyrolysis  area  is  calculated  as  a function  of  time  for  the  three 
different  burning  rates.  These  burning  rates  are  representative  of  unenhanced 
wall  burning  rates  which  would  be  applicable  before  significant  thermal 
feedback  occurs.  The  results  are  shown  in  Figure  4.  The  model  results  with 
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input  burning  rates  of  6 and  8 g/m^s  bound  the  measured  pyrolysis  area  while  7 
g/m^s  compares  well  up  to  700  s.  The  calculated  areas  start  to  level  off 
after  the  pyrolysis  front  reaches  the  ends  of  the  wall  (laterally)  at  the 
wall-ceiling  interface.  The  experimental  results  show  a rapid  increase  in 
pyrolysis  area  due  to  a fast  downward  spread.  This  is  due  to  the  thermal 
feedback  of  the  hot  layer  increasing  the  surface  temperature  of  the  wall,  thus 
increasing  downward  spread.  This  feedback  is  not  explicitly  accounted  for  in 
the  model,  therefore,  we  do  not  expect  to  match  the  experimental  results  at 
long  times.  Also  we  have  used  burning  rates  for  PMMA  associated  with  no 
external  irradiance.  This  is  in  contrast  to  our  use  of  Cone  Calorimeter  data 
at  irradiance  levels  of  25  kW/m^  and  higher  for  corner  fire  tests. 

3.3  Textile  Wall  Covering  Room  Tests 

A series  of  textile  wall  fire  tests  was  performed  by  the  U.C.,  Berkeley  for 
the  American  Textile  Manufactures  Institute  [14].  In  these  tests,  0.31  m (1 
ft)  and  0.62  m (2  ft)  wide  strips  of  different  textile  wall  covering  materials 
were  applied  to  a corner  section  and  the  wall  portion  of  the  wall-ceiling 
interface  of  a room  with  roughly  the  same  dimensions  of  the  room  used  in  the 
tests  above.  A 0.31  m by  0.31  m sand  burner  was  placed  in  the  material  lined 
corner  approximately  10  cm  from  the  walls.  A small  40  kW  fire  was  allowed  to 
burn  in  the  corner  for  5 minutes,  then  the  burner  output  was  increased  to  150 
kW.  The  position  of  the  burner  and  the  size  of  the  small  40  kW  fire  was  such 
that  in  most  cases  it  did  not  cause  the  textile  material  to  ignite,  but 
shortly  after  the  burner  output  was  increased  to  150  kW,  all  wall  lining 
materials  ignited.  The  rate  of  heat  release  was  obtained  during  each  test. 
The  data  for  input  into  the  model  is  given  by  Harkleroad  [15].  Here  lateral 
spread  was  not  included  in  the  calculations  since  lateral  spread  is  not 
expected  with  the  0.31  m wide  strips  anyway.  It  appears  from  the  experimental 
results  that  the  40  kW  fire  only  served  to  preheat  the  room,  therefore  we 
picked  an  initial  surface  temperature  of  100® C based  on  the  maximum  gas 
temperatures  measured  during  this  preheating  time.  In  the  experiments,  there 
appears  to  be  very  little  ignition  delay  of  the  materials  after  the  burner 
output  is  increased  to  150  kW.  This  is  due  to  the  preheating  by  the  small  40 
kW  fire,  and  the  high  heat  flux  from  the  150  kW  fire.  Table  3 shows  the 
experimental  and  calculated  results.  By  examining  the  0.31  m and  0.62  m wide 
strip  results,  it  appears  in  some  cases  lateral  spread  is  significant.  The 
model  calculations  were  performed  for  the  0.31  m wide  strip  cases  using  Cone 
Calorimeter  data  at  external  irradiances  of  30  and  50  kW/m^ . Ignition  delay 
times  for  each  of  the  materials  were  calculated  from  Eq.  (5).  For  this 
calculation,  we  assumed  T^  = 100® C and  = 60  kW/m^  which  represents  the  high 
heat  flux  produced  by  the  ignitor  [16] , and  added  this  calculated  ignition 
time  to  the  model  calculated  time  to  reach  the  peak  rate  of  heat  release.  The 
peak  heat  release  rate  and  the  time  interval  from  the  start  of  the  150  kW  fire 
for  the  0.31  m strips  compare  favorably  with  the  calculations.  After  the  peak 
heat  release  was  obtained  a rapid  decay  was  observed  experimentally  for  most 
cases,  this  was  predicted  well  by  the  model. 
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4.  CONCLDDING  REMARKS 


We  have  presented  evidence  that  an  engineering  approach  utilizing  "fire 
property"  data  from  test  methods  can  satisfactorily  predict  full  scale  energy 
release  of  room  lining  fires.  The  fire  properties  used  are  not  necessarily 
fundamental  engineering  data,  but  the  analysis  does  not  preclude  the  use  of 
such  properties.  Consequently,  the  choice  of  the  fire  properties  used  is 
somewhat  arbitrary  until  the  fire  heat  transfer  conditions  can  be  related  to 
the  test  method  irradiance  levels.  Also  the  model  used  to  make  the  predic- 
tions, although  yielding  reasonably  satisfactory  results,  should  be  refined 
and  improved.  Nevertheless,  the  framework  presented  should  serve  as  a basis 
for  these  improvements  and  might  be  used  now  to  more  rationally  assess  the 
fire  growth  hazard  of  lining  materials. 
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Nomenc latur e 


a 

b 

c 

kpc 

kf 

q 

Q 

t 

t 

X 

y 

r 

n 


area 

dimensionless  parameter,  Eq.  (9) 
dimensionless  parameter,  Eq.  (10) 
dimensionless  parameter,  Eq.  (10) 

thermal  inertia  (thermal  conductivity  x density  x specific  heat) 

constant,  0.01  m^/kW,Eq.  (6) 

heat  transfer 

energy 

time 

temperature 
opposed  flow 
wind-aided  co 
dimensionless 
dimensionless 
dimensionless 
dimensionless 
dimensionless 
dimensionless 
flame  heating 


coordinate  direction 
ordinate  direction 
parameter,  Eq.  (14) 
time , t/tjj 
pyrolysis  position, 
burnout  position,  Xj. 
pyrolysis  position, 
burnout  position,  y^/yg 
parameter,  Eq.  (12) 


> 

/x 


A. 


Subscripts : 


b burnout 

f flame 

ig  ignition 

o initial 

p pyrolysis 


Superscripts : 

(" ) average 
( • ) per  unit  time 
( ) ' per  unit  width 
( ) ” per  unit  area 


11 


5. 


REFERENCES 


1.  Standard  Test  Method  for  Heat  and  Visible  Smoke  Release  Rates  for 
Materials  and  Products  Using  an  Oxygen  Consumption  Calorimeter,  (ASTM  E 
1354),  Amer.  Soc.  for  Testing  and  Materials,  Philadelphia. 

2.  Quintiere,  J.G.  and  Harkleroad,  M. , "New  Concepts  for  Measuring  Flame 
Spread  Properties",  pp.  239-267  in  Fire  Safety  Science  and  Engineering 
(ASTM  STP  882),  Amer. Soc.  for  Testing  and  Materials,  Philadelphia,  PA, 
1985. 

3.  Wickstrom,  U. , Goransson,  U. , "Prediction  of  Heat  Release  Rates  of 
Surface  Materials  in  Large-Scale  Fire  Tests  Based  on  Cone  Calorimeter 
Results",  ASTM  Journal  of  Testing  and  Evaluation,  Vol.  15,  No.  6,  1987. 

4.  Ostman,  B. , and  Nussbaum,  R. , "Correlation  Between  Small-Scale  Rate  of 
Heat  Release  and  Full-Scale  Room  Flashover  for  Surface  Linings",  pp.  823- 
832  in  Fire  Safety  Science  Proc.  of  the  Second  International  Symposium 
(1988),  Hemisphere  Publishing,  New  York  (1989). 

5.  Karlsson,  B. , "Room  Fires  and  Combustible  Linings",  Dept,  of  Fire  Safety 
Engineering,  Lund  Institute  of  Science  and  Technology,  Report,  Lund,  1989. 

6.  Standard  Method  for  Determining  Material  Ignition  and  Flame  Spread 
Properties,  (ASTM  E 1321),  Amer.  Soc.  for  Testing  and  Materials, 
Philadelphia. 

7.  Saito,  F. , Williams,  F.A. , Wichman,  I.S.  and  Quintiere,  J.G.,  "Upward 
Turbulent  Flame  Spread  on  Wood  Under  External  Radiation" , Journal  of  Heat 
Transfer . pp.  438-445,  Vol.  Ill,  May  1989. 

8.  Quintiere,  J.G.,  "The  Application  of  Flame  Spread  Theory  to  Predict 
Material  Performance",  Journal  of  Research  of  the  Nat.  Bur,  of  Stds..  pp . 
61-70,  Vol.  93,  No.  1,  Jan. -Feb.  1988. 

9.  Quintiere,  J.G.,  Harkleroad,  M.  and  Hasemi,  Y. , Combustion  Science 
Technology . Vol.  48,  1986. 

10.  Sundstrom,  B. , "Full  Scale  Fire  Testing  of  Surface  Materials",  Technical 
Report  SP-RAPP  1986:45,  Swedish  National  Testing  Institute,  Boras  1986. 

11.  Quintiere,  J.,  Harkleroad,  M. , "New  Concepts  for  Measuring  Flame  Spread 
Properties",  National  Bureau  of  Standards  (U.S.)  NBSIR  84-2943,  1984. 

12.  Tsantaridis , L.  and  Ostman,  B. , "Smoke,  Gas  and  Heat  Release  Data  for 
Building  Products  in  the  Cone  Calorimeter",  TrateknikCentrum,  Report  I 
8903013,  1989. 

13.  Lee,  B. , Private  Communication,  Memo  to  J.  Quintiere,  Feb.  18,  1987, 
National  Institute  of  Standards  and  Technology,  Gaithersburg,  MD. 


12 


14.  Fisher,  F. , MacCracken,  B.,  and  Williamson,  R. , "Room  Fire  Experiments  of 
Textile  Wall  Coverings",  Fire  Research  Laboratory,  University  of 
California,  Berkeley,  California,  1986. 

15.  Harkleroad,  M. , "Fire  Properties  Database  for  Textile  Wall  Coverings", 
National  Institute  of  Standards  and  Technology  (U.S.),  NISTIR  89  4065, 
1989. 

16.  Mowrer,  F. , Private  Communications,  Univ.  of  Maryland,  July  1990. 


13 


Table  1 


Flame  Spread  and  Heat  Release  Properties  of  Swedish  Fire  Test  Materials 


T. 

1 s 

kpc 

$ 

, min 

Q" 

Q" 

Material 

(“C) 

(kW/m^K)^  s 

(kW^/m^  ) 

(“O 

(kW) 

(25 

(s) 
kW/m^ ) 

(kW)  (s) 

(50  kW/m^) 

Particle  Board 

405 

0.626 

8 

180 

140 

>500 

200 

>500 

Insulating 

Fiberboard 

381 

0.229 

14 

90 

120 

>500 

150 

>300 

Medium  Density 
Fiberboard 

361 

0.732 

11 

80 

140 

>500 

170 

>600 

Wood  Panel 
(Spruce) 

389 

0.569 

24 

155 

140 

>200 

160 

>400 

Melamine  Covered 
Particle  Board 

483 

0.804 

<1 

435 

- 

- 

115 

>400 

Paper  Covered 
Gypsum  Board 

388 

0.593 

0.5 

300 

100 

40 

200 

20 

PVC  Covered 
Gypsum  Board 

410 

0.208 

25 

300 

105 

20 

160 

20 

Textile  Covered 
Gypsum  Board 

406 

0.570 

9 

270 

200 

30 

300 

20 

Textile  Covered 
Mineral  Wool 

391 

0.183 

6 

174 

250 

30 

375 

20 

Paper  Covered 
Particle  Board 

426 

0.680 

13 

250 

140 

>500 

150 

>500 

Polyurethane 

Foam  (rigid) 

393 

0.031 

3 

105 

110 

60 

180 

60 

Expanded 

Polystyrene 

482 

0.464 

- 

- 

- 

- 

300 

100 

Gypsum  Board 

469 

0.515 

14 

380 

- 

. 

100 

20 

Data  were  not  taken 
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Table  2 


Time  to  Achieve  a 1 MW  Fire  for  the  Swedish  Room  Tests 


Material 


Experimental  Time 
(s) 


Model  Calculations 


Irradiance 
25  kW/m2 
(T3  = 25-C) 
(s) 


Irradiance 
50  kW/m^ 
(T3  = 25°C) 
(s) 


Irradiance 
25  kW/m^ 
(T3  = 80"C) 
(s) 


Particle  Board 

157 

198 

143 

145 

Insulating 

Fiberboard 

59 

77 

58 

55 

Medium  Density 
Fiberboard 

131 

180 

148 

125 

Wood  Panel 
(Spruce) 

131 

165 

143 

117 

Melamine  Covered 
Particle  Board 

465 

- 

402 

- 

Paper  Covered 

Gypsum  Board 

640 

632 

616 

641 

PVC  Covered 

Gypsum  Board 

611 

619 

606 

622 

Textile  Covered 
G3rpsum  Board 

639 

615 

613 

615 

Textile  Covered 
Mineral  Wool 

43 

33 

28 

24 

Paper  Covered 
Particle  Board 

143 

237 

220 

177 

Rigid  Polyurethane 
Foam 

6 

11 

6 

7 

Expanded 

Polystyrene 

115 

- 

122 

- 

Gypsum  Board 

ic 

- 

■k 

- 

Data  were  not  taken 
* Did  not  Reach  1 MW 
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Table  3 


Textile  Wall  Coverings  Room  Fire  Tests 


Full 

Scale 

Screening  Tests 

Model 

Calculations 

0, 

.31  m 

Strips 

0.62m 

Strips 

Irradiance 

30  kW/m^ 

Irradiance 

50  kW/m^ 

Qp 

Qp 

Qp 

Qp 

(kW) 

(s) 

(kW) 

(S) 

(kW) 

(s) 

(kW) 

(s) 

MATERIAL 

(H)  85%  wool 

15%  cotton 

46 

30 

160 

40 

- 

146 

46 

(C)  55%  cotton 
45%  rayon 

62 

30 

119 

60 

137 

40 

139 

37 

(G)  100% 
polyester 

83 

30 

- 

64 

39 

56 

44 

(B)  100% 
polyester 

207 

45 

298 

60 

121 

46 

270 

46 

(Q)100% 

polyester 

207 

40 

480 

40 

145 

50 

293 

55 

(Qfr)  100% 
polyester 

310 

40 

- 

157 

43 

292 

59 

(R)  100% 
nylon 

587 

70 

590 

70 

46 

46 

416 

51 

(AA)  70%  acrylic 
30%  wool 

684 

30 

- 

725  109 

744 

106 

(PPPF) 

- 

337 

50 

271 

45 

450 

48 

polypropylene 


- Data  were  not  taken  Qp  = peak  energy  release  rate 

t*  = time  interval  from  start  of  150  kW  burner  to  peak  energy  release  rate 
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Figure  1.  Energy  release  property  data. 
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Figure  2.  Mode!  floor,  ceiling  and  wall  configurations 
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FULL  SCALE  ROOM  CORNER  TESTS 


(M>i)  3SV3"I3a  1V3H  30  31Va 


Figure  3.  Energy  release  predictions  for  the  Boras  data 
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Figure  4.  Growth  of  the  pyrolysis  area  on  a PMMA  wall 
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Appendix  A.  Full  Scale  Room-Corner  Tests 


The  experimental  results  and  model  predictions  for  the  twelve  materials 
examined  at  the  Swedish  National  Testing  Institute  are  shown  in  figures  A-1  to 
A-12.  Model  calculations  with  T^  = 25*C  and  SO^C  are  both  presented.  Data 
From  the  Cone  Calorimeter  at  an  irradiance  of  25  kW/m^  was  used  for  these 
calculations  with  the  exception  of  Melamine- faced  particle  board  and  Polys- 
tyrene where  Cone  Calorimeter  data  at  an  irradiance  of  50  kW/m^  were  used. 
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RATE  OF  HEAT  RELEASE  (kW)  RATE  OF  HEAT  RELEASE  (kW) 


0 300  600  900 

TIME  (S) 

Figure  A-1<,  Rate  of  heat  release  vs.  time  for  particle  board 
lined  room 


0 300  600  900 

TIME  (S) 

Figure  A-2.  Rate  of  heat  release  vs.  time  for  insulating 
fiberboard  lined  room  fire 
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RATE  OF  HEAT  RELEASE  (kW)  RATE  OF  HEAT  RELEASE  (kW) 


0 300  600  900 

TIME  (S) 


Figure  A-3.  Rate  of  heat  release  vs.  time  for  medium  density 
fiberboard  lined  room  fire. 


0 300  600  900 

TIME  (S) 


Figure  A-4.  Rate  of  heat  release  vs.  time  for  wood  paneling 
lined  room  fire. 
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RATE  OF  HEAT  RELEASE  (kW)  RATE  OF  HEAT  RELEASE  (kW) 


1500 


1 000 


500 


300 


TIME  (s) 


600 


900 


Figure  A-5.  Rate  of  heat  release  vs.  time  for  melamine  faced 
particle  board  lined  room  fire. 


0 300  600  900 

TIME  (S) 

Figure  A-6.  Rate  of  heat  release  vs.  time  for  paper  covered 
gypsum  board  room  fire. 
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RATE  OF  HEAT  RELEASE  (kW)  RATE  OF  HEAT  RELEASE  (kW) 


0 300  600  900 

TIME  (S) 

Figure  A-7,  Rate  of  heat  release  vs.  time  for  pvc  covered 
gypsum  board  room  fire. 


TIME  (S) 


Figure  A-8.  Rate  of  heat  release  vs.  time  for  textile  wallcovering 
on  gypsum  board  lined  room  fire. 
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Figure  A-9-  Rate  of  heat  release  vs.  time  for  textile  wallcovering 
on  mineral  wool  lined  room  fire. 
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Figure  A-10.  Rate  of  heat  release  vs.  time  for  paper  covered 
particle  board  lined  room  fire. 
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Rgure  A-11.  Rate  of  heat  release  vs.  time  for  rigid 
polyurethane  foam  lined  room  fire. 
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Figure  A-12.  Rate  of  heat  release  vs.  time  for  expanded 
polystyrene  lined  room  fire. 
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Appendix  B.  PMMA  Room  Wall  Fire  Test 


Figure  B-1  shows  surface  temperature  measurements  for  a PMMA  room  wall  fire 
performed  at  NIST  [12].  The  (x)  symbols  denote  the  locations  of  the  0.25  mm 
chromel-alumel  thermocouples,  the  (o)  symbols  denote  the  locations  of  heat 
flux  gages.  The  numbers  indicate  the  time  in  seconds  from  the  start  of  the 
test  to  the  time  when  the  thermocouple  reads  a surface  temperature  of  378® C. 
This  temperature  is  used  as  the  critical  surface  temperature  at  ignition.  The 
underlined  numbers  indicate  the  time  in  seconds  for  the  heat  flux  gage  to 
register  a heat  flux  of  15  kW/m^  . Table  B-1  is  an  estimate  of  pyrolysis  area 
as  a function  of  time  for  this  experiment.  These  values  are  plotted  in  Figure 
4 as  the  experimental  values . 


28 


Table  B-1 


Pyrolysis  area  estimated  from  thermocouple  data 


Time 

(s) 

Pyrolysis  Area 
(m2) 

245 

345 

435 

520 

580 

675 

700 

750 

800 

0.19 

0.28 

0.40 

0.47 

0.56 

1.39 

1.95 

3.62 

4.18 
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Figure  B-1.  Time  to  reach  ignition  temperature  at  various 
positions  for  flame  spread  on  a PMMA  wall  section  in  an 
enclosure. 
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